Spin accumulation in ferromagnetic single-electron transistors 
in the cotunneling regime 
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We propose a new method of direct detection of spin accumulation, which overcomes problems of 
previous measurement schemes. A spin dependent current in a single-electron transistor with ferro- 
magnetic electrodes leads to spin accumulation on the metallic island. The resulting spin-splitting 
of the electrochemical potentials of the island, because of an additional shift by the charging energy, 
can be detected from the spacing between two resonances in the current-voltage characteristics. The 
results were obtained in the framework of a real-time diagrammatic approach which allows to study 
higher order (co-)tunneling processes in the strong nonequlibrium situation. 
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The discovery of materials and devices with spin de- 
pendent electronic transport properties opens new per- 
spectives for research and applications. As examples 
we mention exchange-coupled metallic magnetic multi- 
layers, heterostructures with ferromagnetic semiconduc- 
tors, hybrid structures based on magnetic metals and 
nonmagnetic semiconductors, as well as perovskite ox- 
ides. The field has developed into a new branch of meso- 
scopic electronics - called magneto- or spin electronics 
E[. The balance of spin injection and spin diffusion in 
spin-polarized heterostructures in general leads to a non- 
equilibrium spin accumulation |^||, which in turn pro- 
duces a difference in the local electrochemical potentials 
for spin-up and -down electrons. The predicted spin po- 
larization was probed in experiments, either by measur- 
ing the resistance of metallic devices |^,^| or in semicon- 
ducting systems || by studying the circular polarized 
photoluminescence . 

In spite of the reported successes in the indirect de- 
tection of spin accumulation, there remain controversies 
about the results || . In these methods one usually mea- 
sures the resistance and for a typical experimental config- 
uration some other effects related to interface scattering, 
anisotropic magnetoresistance, etc, can occur. Even the 
question whether in principle it is possible to observe the 
spin splitting of the electrochemical potential by spectro- 
scopic methods (analogous to the tunneling spectroscopy 
for the superconducting gap) remains open. The reason 
is the fact that the spin accumulation is a nonequilib- 
rium effect, and the splitting of electrochemical potential 
is always smaller than the applied bias (in contrast to 
the superconducting energy gap which does not depend 
on the applied voltage). 

In this article, we propose a method of direct detect- 



ing the nonequilibrium spin splitting of the electrochemi- 
cal potentials in ferromagnetic single-electron transistors 
(FM SETs), which is free of the difficulties met in the 
methods used previously. We show that the splitting can 
be read from the spacing between two resonances in the 
current-voltage characteristics. The resonances are sepa- 
rated due to the spin dependence of the electrochemical 
potentials, but in addition they are shifted by the charg- 
ing energy due to Coulomb blockade effects. 

Spin-dependent transport in ferromagnetic double- 
barrier tunnel junctions has been studied before, both ex- 
perimentally M and theoretically in the sequential ^ |o) , 
cotunneling |Tl]Jl^| and strong tunneling jl3| regimes. 
The previous work covering the cotunneling regime is 
based on the approach developed by Averin and Nazarov 
Jhij for nonmagnetic SETs, which, however, is valid only 
away from the resonance. On the other hand, exactly 
the resonance is crucial for the effects considered in this 
paper. Therefore, we extend the previous descriptions 
by using the diagrammatic real-time technique, devel- 
oped for normal-metal SETs |l5|Jl^ ], and extend it to 
FM SETs with spin-dependent tunneling. We take into 
account single-barrier cotunneling processes p8[ , vertex 
correction and propagator renormalization. 

The system under study consists of a small nonmag- 
netic metallic grain (island) connected to two ferromag- 
netic leads (electrodes) via tunnel barriers |p| Jlfj| , p^| . Its 
Hamiltonian takes the form 
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states of the leads and island are described by the wave- 
vectors k and q, transverse channel index u, and spin a. 
The Coulomb interaction on the island is accounted for 
by f/ch = Ec(n — n x ) 2 , where Ec = e 2 /2C is the scale 
for the charging energy, and C = Cl + Cr + C g is the 
total capacitance of the island, which is the sum of the 
left and right junction capacitances and the gate capaci- 
tance. The 'external charge' en x = ClVl + CuVR + C g V g 
accounts for the effect of the applied voltages, Vl and Vr 
in the left and right electrodes and the gate voltage V g . 
The last part of the Hamiltonian, £/t, describes tunnel- 
ing processes and may be written as 



r=L,R kq no 



(2) 



The phase operator (p is the conjugate to the charge eh 
on the island, and the operator e ±l ^ describes changes of 
the island charge by ±e. 

When writing Eq. (0) we assumed that the electron 
spin and transverse channel index are conserved dur- 
ing tunneling. We further will assume for simplicity 
that the transfer matrix elements T^ v are only de- 
pendent on the junction r and the spin orientation a, 
Tf.° — T*. They can be related to the spin depen- 
dent tunneling resistance of the barriers via the relation 
1/Rra = {2TTe 2 /fi)N\T 1 a \ 2 D la D u where D ra is the spin- 
dependent density of electron states at the Fermi level in 
the electrode r, D\ the density of states per spin in the 
island, and N the number of transverse channels. 

We first ignore intrinsic spin-flip processes in the is- 
land, since usually the spin flip time r s f is much longer 
than the time scale for the variation of the charge on the 
island and the energy relaxation time. Thus, we treat the 
two spin components of the electrons on the island as two 
independent electron reservoirs. They are characterized 
by Fermi distributions fi a (E) = l/{exp[/3(E — fii a )] + l}, 
where are the corresponding spin dependent electro- 
chemical potentials which are different in a general case, 
MiT Mi| ■ Without loss of generality we can choose 
the reference energy such that they satisfy //if = — //ij.. 
Spin current conservation, as described later, fixes these 
values. 

The grand-canonical density matrix of the system de- 
pends on the electrochemical potentials of the two spin 
components. Following the standard procedure [ p^|Jl6| 
we expand it in ifr(i) and perform the trace over the 
reservoir degrees of freedom using Wick's theorem. The 
electric current is given in the lowest order perturbation 



expansion by = J2 a I } 
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Here A n = E c ^{n + 1) — E c h(n) , and af a {e) are the for- 
ward and backward propagators on the Keldysh contour 
in Fourier space, 



a±(e)=±a r ° 



exp[±/3(e - A//„)] - 1 
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The dimensionless conductance of the junction r for spin 
a is a CT = h/(4:iT 2 e 2 R la ), e is the energy of the tunneling 
electron, and A/i r(T = fj, T — \i\ a . Apart from this, a (e) = 
Y^rcr a t<T( e )i a ( € ) — a+ i e ) + a ~( e )> an d the probabilities 



Pn^ obey the equation pn'a + (A n ) — Pn +1 a (A„) = 

witn E n P™ = l - 

The dominant second order (cotunneling) contribution 
to the electric current can be divided into three parts, 

' (2) = EL 4^. with jf> = £,i«= r £ c r(2) 

These terms are given by the following equations: 
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where #„(u>) = l/(w - A n + i0+) - 1/(cj - A n _! + i0 + ) . 



The terms and describe renormalization of the 
tunneling conductance and energy gap, respectively, and 
become important at resonance. 

We have numerically evaluated the current for par- 
allel (P) and antiparallel (AP) alignment of the elec- 
trode magnetizations. We assumed symmetric junctions 
such that for parallel alignment Rr-\ = i?L| = &Rq & n d 

R*i = Ru = (i - P)/0- + P)RRh where Rq = 

and we assume a = 5. The parameter P denotes the 
spin polarization of the electrodes - P = 0.23 and 0.40 
for Ni and Fe, respectively. For the antiparallel align- 
ment one then finds = i?Li = clRq, Rri — Rli = 
(1-P)/(1 + P)R RV _ 

The spin accumulation on the island (or equivalcntly 
spin splitting of the electrochemical potential) is deter- 
mined from the condition of spin current conservation 



r (2) 



(lg>+lW)-en <r D l /T B[ = 0. 
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Here the spin flip processes in the island are taken into 
account, and characterized by the relaxation time r B f. 

Equation (|J) includes the components, which describe 
single-barrier cotunneling processes which transport spin 
but do not transport charge. They cancel after sum- 
mation over the spin orientation and do not enter the 
current formula. On the other hand, they enter Eq. (||) 
and, by modifying the magnetic state of the island, have 
an indirect influence on the transport current. 

In Fig. [I] we show the differential conductance for Ni 
and Fe electrodes in both parallel and antiparallel mag- 
netic configurations, calculated for T/Eq = 0.02 in the 
absence of spin-flip processes. The corresponding value 
of tunnel magnetoresistance, TMR = (Rap — Rp)/Rp 
where Rp and Rap are total resistances for parallel 
and antiparallel configuration respectively, is shown in 
Fig. ||(a). In the AP configuration the splitting of the 
conductance peak, resulting from the spin splitting of the 
electrochemical potential of the island, is well resolved 
(see Fig. 0(b)). Different splitting for Fe and Ni is due 
to different magnetic polarizations P of these metals. 

The magnitude of the associated spin accumulation is 
shown in Fig. §(b). It is interesting to note that no spin 
accumulation, and consequently no splitting of the reso- 
nance, occurs in the P configuration (see Fig. [j](a)). 

The cotunneling component of the current shows reso- 
nances when the bias voltage V approaches the - effective 
- electrochemical potential of either spin component. The 
pure spin splitting of the electrochemical potentials is al- 
ways smaller then the bias voltage. However, Coulomb 
blockade effects effectively shift the electrochemical po- 
tentials further by the charging energy (as illustrated in 
the diagram of Fig. ||). Due to this additional shift the 
resonances in the current should be resolvable. A direct 
observation of the resonances in the current may still 
be difficult because even at resonance cotunneling yields 
only a small fraction of the total current. The resonances 
arc easier to detect in the differential conductance, which 
is shown in Fig. |l|. For the interpretation of these results 
it should be noted that a resonance in the current corre- 
sponds to a maximum negative slope of the conductance. 

Such an effect may also arise in the sequential tunnel- 
ing limit p0|| (high resistance junctions), but in that case 
it is not so evident and pronounced as in the cotunneling 
limit (see Fig. ^). The role of the cotunneling current 
is twofold: (i) the cotunneling current has a sharp maxi- 
mum at resonance (described above) , which can be easily 
detected on the conductance characteristics (see compar- 
ison of cotunneling and sequential limit in Fig. ^), and (ii) 
in the Coulomb blockade regime sequential tunneling is 
suppressed, so only the cotunneling processes contribute 
to spin accumulation (as can be concluded from Fig. ^|(b) 
for eV/2E c < I). 

In Fig. | we show the differential conductance versus 
gate V g and transport V voltages in a gray-scale plot for 
both P and AP magnetic configurations. The splitting of 



the resonance (shown in Fig. |l|(b) for V g — 0) is clearly 
visible in Fig. ||(b). From this figure one can read out 
the dependence of the electrochemical potential splitting 
on the transport voltage V. Figure |i|(b) also shows that 
the splitting of the conductance peaks can be observed 
when the gate voltage V g is varied and the bias voltage 
V is fixed. 

Generally, spin-flip relaxation processes suppress spin 
accumulation. Figure |5| shows the values of the volt- 
age corresponding to the upper and lower resonances 
from Fig. |l](b) (maximum negative slopes of the con- 
ductance curves) as a function of the intrinsic spin re- 
laxation time r s f on the island. Generally, the spin- 
relaxation time r s f depends on sample properties. For 
instance, in Al one finds r s f ~ 10~ 8 s. For small islands, 
D\ ~ 10 3 /eU, and for the other parameters as in Fig. [I], 
one finds r s fi?Q//i(i?L + Rr)Dj ~ I0 3 . This means that 
the spin relaxation process is relatively slow and the value 
of the spin accumulation is close to its maximum value 
for Tgf — * oo. It is known that spin accumulation can 
be suppressed by a small perpendicular magnetic field 
('Hanle effect') M. In principle, this effect can be mea- 
sured in the system under consideration and may be used 
to determine the spin relaxation time as well. 

There are several experimental works |7j on ferromag- 
netic SETs both in the sequential and cotunneling regime 
with all electrodes (including island) ferromagnetic. The 
corresponding data, as well as the data obtained on fer- 
romagnetic granular systems suggest that spin accumu- 
lation may occur in such systems, but there was no clear 
experimental evidence on this. But the range of the pa- 
rameters indicates a possibility of a direct detection by 
the method proposed in this paper. 

In conclusion, we have analyzed electron tunneling 
in ferromagnetic SETs in the cotunneling regime. We 
showed that the spin accumulation on the island can be 
observed directly in the current-voltage characteristics as 
a splitting of the resonance. The distance between two 
resonances arises due to the spin-splitting of the electro- 
chemical potential due to spin accumulation. This effect 
should be observable experimentally, and could give clear 
evidence of the spin accumulation, overcoming the objec- 
tions to the previous measurements . 
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FIG. 1. Differential conductance vs. bias voltage V in the 
(a) parallel and (b) antiparallel configurations, calculated for 
T/Ec = 0.02 and symmetric junctions with i?R| = 5-Rq. The 
spin polarization is P = 0.23 and 0.40 for Ni and Fe electrodes, 
respectively. Dotted line is for the sequential tunneling limit 
for Fe electrodes. 
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FIG. 2. (a) Tunnel magnetoresistance TMR and (b) spin 
accumulation for Ni and Fe electrodes as functions of bias 
voltage V. The other parameters are the same as in Fig. 
Dotted line is for the sequential tunneling limit for Fe elec- 
trodes. 
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FIG. 3. Energy diagrams for a symmetric magnetic dou- 
ble tunnel junction with a normal metallic island for the an- 
tiparallel configuration (a) without and (b) with the Coulomb 
blockade. Here //l, Mr are the electrochemical potentials for 
the left and right electrodes and fj,f , f±i are the electrochemical 
potentials for spin-up and spin-down electrons on the island. 
N = 1 , denotes the island charge. 
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FIG. 4. The differential conductance vs. gate V g and trans- 
port V voltages in a gray-scale representation in the (a) par- 
allel and (b) antiparallel configurations calculated for Fe elec- 
trodes and the other parameters as in Fig. nj. 
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FIG. 5. Dependence of the upper and lower resonance volt- 
ages on the intrinsic spin relaxation time r 8 f, calculated for Fe 
electrodes and the other parameters as in Fig. |l]. The dotted 
line denotes a regime of insufficiently resolved resonances. 
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